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ABSTRACT 

We describe a powerful signal processing method, the continuous wavelet transform, 
and use it to analyze radial structure in Cassini ISS images of Saturn's rings. Wavelet 
analysis locally separates signal components in frequency space, causing many struc- 
tures to become evident that are difficult to observe with the naked eye. Density waves, 
generated at resonances with saturnian satellites orbiting outside (or within) the rings, 
are particularly amenable to such analysis. We identify a number of previously unob- 
served weak waves, and demonstrate the wavelet transform's ability to isolate multiple 
waves superimposed on top of one another. We also present two wave-like structures 
that we are unable to conclusively identify. 

In a multi-step semi-automated process, we recover four parameters from clearly 
observed weak spiral density waves: the local ring surface density, the local ring viscosity, 
the precise resonance location (useful for pointing images, and potentially for refining 
saturnian astrometry), and the wave amplitude (potentially providing new constraints 
upon the masses of the perturbing moons). Our derived surface densities have less 
scatter than previous measurements that were derived from stronger non- linear waves, 
and suggest a gentle linear increase in surface density from the inner to the mid-A Ring. 
We show that ring viscosity consistently increases from the Cassini Division outward to 
the Encke Gap. Meaningful upper limits on ring thickness can be placed on the Cassini 
Division (3.0 m at r ~ 118,800 km, 4.5 m at r ~ 120, 700 km) and the inner A Ring 
(10 to 15 m for r < 127, 000 km). 
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Introduction 



The radial structure of Saturn's rings has long been the subject of inquiry by astronomers, 
dating from G. D. Cassini's discovery in 1675 of the division that bears his name, and W. Herschel's 
realization in 1791 that the Cassini Division appears identica l on both the nort hern and southern 
ring faces, definitively identifying it as a gap (lAlexandeiill962l : IVan Heldenlll984l ). Many additional 
ringlets and gaps came to light during succeeding centuries, but more subtle variations of density 
within a ring awaited discovery by spacecraft. 

In preparation for spacecraft encounters, an analogy with galactic dynamics - pa rticularly the 
Lindblad resonances and spiral density waves that cause galaxies to have spira l a rms |Lin and Shu 
1964 ) - was first applied to Saturn's rings by Goldreich and Tremaine ( 1978al lbl. Il980 ). who noted 
that radial structure could arise due to the perturbations of e xterior satellites. The general th eory 
of resonant in teractions within planetary rings was reviewed bv I Goldreich and Tremaind (|l982l ) and 
byH |l983). 



The Voyager encounters with Saturn in 1980 and 1981 produced an u nprecedented data 



set on the rings' radial structure, including evidence for spiral density waves (jCuzzi et al. 



1981 



Holberg et a/.lll982l : llIolberall982l ). A number of authors quantitatively analyzed the radial traces 
of density waves; most of these focused on deriving the ring's local backg round surface density , 
which is directly obtained from the wavelength dispersion of density waves. Esposito et al. (jl983l ) 
carried out the first systematic survey for wave traces, firmly establ ish ing the connect ion with the 



known locations of Lindblad resonances (jLissauer and Cuzzil Il982l ). lEsposito et al\ additionally 
attempte d to fit the waves' dam ping length and thus derive ring viscosity. The analysis of wave 
traces by [Nicholson et al\ (|l99Cl ) concentrated on pinpointing the locations of resonant wave gener- 
ation, aiding thei r purpose of precisely navigating the occultation scans and determining Saturn's 
pole orientation. Rosen et al. ( 1991al lb[) wrote extensively on methods of density wave analysis, 
and addition a lly fit wave amplitudes in order to derive masses for the perturbing moons. Finally, 
Spilker et al\ (|2004l ) undertook a comprehensive survey of spiral density waves in the Voyager data 
set. 

A distinct form of radial structure in the rings are the wakes raised by the gravity of a pass- 
ing nearby moon (here called "moonlet wakes" to distinguish them from the "sel f-gravity wakes" 
menti oned briefly in Section 16. 2p . These were described and explained in detail by IShowalter et al. 



(j 19861 ). who used that under standing to pinpoint and then obser ve the moon Pan embedded in the 



Encke Gap (jShowalteii Il99ll ) . Further analysis of Pan wakes by iHorn et al\ (|l996l ) was among the 
first to use a form of loc alized Fourier analysis to examine the rings' radial structure (though see 
also lEsposito et al.lll983l). The sam e technique was applied to the unexplaine d radial structure i n 



the B Ring (jHorn and Cuzzilll996l ) and to spiral density waves in the A Ring (jSpilker et a/.ll2004l ). 



The Imaging Science Subsystem (ISS) (jPorco et a/.ll2004l ). on board the Cassini spacecraft, has 
vastly improved upon the clarity of Voyager images of the rings. Cassini's highest-quality imaging 
to date of Saturn's rings was accomplished during the spacecraft's insertion into Saturn orbit (SOI) 
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on 1 July 2004. The companion paper by iPorco et al\ (|2007l . hereafter iPaper 11 ) describes the 
cahbration, image processing, and presentation of these data, as well as of images collected during 
a lower-resolution but contiguous radial sequence across the sunlit side of the rings on 20 May and 
21 May 2005, resulting in a series of scans of brightness with orbital radius. Also in iPaper 11 is a 
catalog of the strongest gravitational resonances and resonance strengths arising from the known 
moons, including Pan and Atlas, both within and exte rior to the rings, derived using the latest 



Cass ini values for moon masses and orbital elements (jJacobson et al. 



2005 



2006 : ISpitale et al 



In this paper, we analyze the 1 July 2004 radial scans usiiig techniques involvin g the wavelet 
transform (jDaubechied Il992l : iFargd Il992l : iTorrence and Compd Il998l : lAddisonl |2002| ) , which func- 
tions as a method of localized Fourier analysis that is optimized to simultaneously d etect signals 



over a wide range of frequencies. Developing the wavelet techniques first presented by lBurns et al 



(j200J), we employ a multi-step semi-automated process to measure wave parameters that can be 
inverted to yield ring and moon properties. In this task, we focus on density waves that are too 
weak to have been detected by Voyager, not only to avoid du plication of previous w ork but also 
because these weaker waves maintain strictly linear dispersion (|Shu et a/.lll983l . Il985l ). causing our 
results to be less susceptible to systematic errors. Additionally, we present several examples of 
wavelet transforms of the SOI scans, both to further illuminate the structure of Saturn's rings and 
to demonstrate the advantages of wavelet analysis. 

A brief paper describing our model of the complex spiral density wa,ves ge ner a,ted at reso nances 
with the co-orbital satellites Janus and Epimetheus (jTiscareno et al\ l2006l ) is iPaper III] in this 
sequence. 

Section [2] reviews the theory of spiral density waves. Section [3] that of wavelet transforms, and 
Section H] describes a new wavelet method for analyzing density waves. Section [5] shows examples 
of Cassini ISS radial scans of the rings, with features explicated by wavelet analysis. Section [6] 
presents the results of our analysis of weak linear density waves. Section [7] presents our conclusions. 



2. Density Waves 



We begin by recalling the mathematical properties of spiral density waves. Such waves occur 
at locations wher e ring particles are in an inner or outer L indblad resonance (ILR or OLR) with a 
perturbing moon ([Goldreich and Tremainelll982l : IShulll984l ). At such locations, the frequency of the 
moon's perturbation forms a simple ratio with the radial (epicyclic) frequency of ring particle orbits, 
amplifying their eccentricity. The coherent excitation of ring particles gives rise to compression and 
rarefaction, initiating a wave that propagates through the disk. 

The density wave generated in a ring at an inner Lindblad resonance (ILR) can be described 
as a density variation Act on the b ackground surface density an. At orbital radius r greater than 



the resonance location r^, we have ( Goldreich and Tremaine 



1982 



Shu 1984: Nicholson et al. 



199C 
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Rosen et a/.lll991al fi] 



Aa{r) = Re < iAie 



-i(j>0 



1 - i^e 



where the dimensionless radial parameter is 



\2TrGao 



-ie/2 



1/2 



drj 



(1) 



rL 



(2) 



and further terms are defined below. 



Assumi ng Saturn's gravity is well described as a p oint mass plus a J2 harmonic, the factor SIl 
is given by (ICuzzi et a/Jll984l : iMarley and Porcdll993l ) 



= 3(m - + J2 



21 

T 



9 



(m — 1) 



02 



(3) 



where the second term is a small correction except for m = 1. 



The Lindblad resonance occurs at r^. The local mean motion is r?/-,, which must be calculate d 
accounting for the higher-order moments of Sat urn's gravity field (jMurrav and DermottI Il999l ). 
Saturn's radius rs = 60,330 km by convention (jKliore et a/.lll98d ). The resonance's azimuthal 
parameter is m, a positive integer equal to the number of spiral arms. The amplitude Al is related 
to the mass of the perturbing satellite, while the damping constant describes the ring's viscous 
response (see Section [2?3l) . 



2.1. Initial Phase 

For a given longitude on the rings. A, the density wave's initial phase is given by 

(pQ = mX — {m + k)\s + kwg, (4) 

where A^, and Ws are the mean longitude and longitude of periapse of the perturbing satellite. 
The integers m and k describe the resonance — a (/i;+l)th-order Lindblad resonance is commonly 
labeled {m+k):{m-l) . A third integer p describes resonances linked to the perturbing satellite's 
inclination, but no such resonances are discussed in this paper. 



'There are differences in notation among these references. Fi rstly, [Sh^ ill 9841') and [Nicholson et al. I lIlQQd . hereafter NCP) 
define ^ to be a factor of \/2 smaller than our Eq. [2] (which follows fRosen et al\\l991s^ ). Secondly, all three papers use diff'erent 
normalizations and absorb them into the amplitude Aj^ . Specifically, 



— ~i o~ — ^O^L 
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2.2. Local Wavenumber and Phase 



The integral in Eq. [T] is a Fresnel integral, which significantly modulates the result near the 
wave's generation point, but oscillates about unity for higher values of ^. Downstream, then, the 
dominant component of Eq. [1] has the form of a sinusoid with constantly decreasing wavelength 
(as well as modul ating amp litude), such that the wavenumber, k = 27r/A, increases linearly with 
distance from ri (|Shulll984l ): 

kDw{r) = ^ ^ (r-r^). (5) 



2TrGaorL 



Since the wavenumber increases linearly with r, the accumulating phase naturally increases 
quadratically. Indeed, making use of J kdr = ^^/2, the instantaneous phase of the sin usoid far fr om 
^ = (i. e., r = r^), as a function of radial location ■^(r), asymptotically approaches (jShulll984l ) 



</>w = <Ao + f/2 + 7r/4. 



It is important to note the term of 7r/4 (or 45°) added to this quantity. This asymptotic limit 
good approximation only for ^ ^ 1 (jNicholson et a/.lll99Cll ). 



(6) 



IS a 



2.3. Ring Viscosity and Vertical Scale Height 



When the damping constant ^ 1 (that is, when damping is inefficient enough to allow at 

least s e veral wave cycles), the local ring kinematic viscosity u can be estimated as (jGoldreich and Tremaine 



1978bl : IShulll984l ) 



9 



1/2 



(7) 



where, for this purpose, the radial frequency is approximately equal to the orbital frequency 
The viscosity v is directly related to the rms random velocity c ( Goldreich and Tremainjl978a ) 



(8) 



where r is the local optical depth. Eq. [8] assumes that ring particle interactions are isolated 
two-particle collisions; it is not valid when the ring p article density is high enough that particle 
size becomes important ( Wisdom and Tremainelll988l ). This assumption is probably valid for the 
A Ring. 



Under the assumption that random velocities are isotropic, the ring's vertical scale height can 
be estimated as if ~ c/J7. However, this assumption is violated in much of the A Ring (Section [62]) • 
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2.4. Summary 

The density wave equation (Eq. [1]) contains five parameters that may vary from wave to wave 
(not counting the azimuthal parameter m, which is fixed by the identity of the resonance), and 
which are sensitive to physical quantities concerning the ring or the perturbing moon. 

1. The background surface density of the ring, ctq, governs the rate at which the wavenumber 
increases with distance from resonance. 

2. The resonance location, ri, fixes the wave against translation in the radial direction. 

3. The initial phase, (f)Q, is fixed by the moon's orbital phase relative to the location at which 
the wave is observed. 

4. The damping parameter, ^^i, governs the location at which the wave's amplitude ceases to 
grow and begins to decay; it is sensitive to the dynamical viscosity of the ring. 

5. The amplitude, Al, governs the overall strength of the perturbation (irrespective of its shape); 
it is sensitive to the mass of the perturbing moon. 



3. Wavelet Theory 



This paper presents the first thorough analysis of periodicities in the radial structure of Saturn's 
rings using the wavelet transform, a powerful method of signal processing that has been applied 
successfully in many other fields. Since this method is complex, we review its properties in this 
section, emphasizing aspects that are relevant to our purpose s. For more detailed information on 
the wavelet transform and its applications, see lAddisonl (120021 ) and references therein. 



3.1. The Morlet Wavelet 

The mother wavelet governs the way in which the wavelet transform transfers spatial infor- 
mation into the frequency domain. The mother wavelet we use is the Morlet wavelet, which is 
optimally suited for identifying oscillatory components of a signal. 

A wavelet transform using the Morlet mother wavelet can be thought of as a localized Fourier 
transform. As with standard Fourier analysis, components of the input signal are separated in the 
frequency domain, but wavelet analysis additionally pinpoints the locations in the signal at which 
a given frequency is important. This is particularly useful when the signal's dominant frequency 
changes with location (as it does in a density wave propagating through a disk). The advantage of 
wavelet analysis over other time-frequency methods (such as the windowed Fourier transform) lies 
in the automatic scaling of the window along with the interrogating waveform, yielding optimal 
simultaneous detection of both high- and low-frequency signal components. 
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The Morlet wavelet is simply a complex sinusoid within a Gaussian envelope: 

il;{t) = TT-^/"^ exp [iujot - . 
The wavelet's central frequency, wq, is selected beforehand by the user (see Section 



(9) 



Wavelet analysis is carried out by translating and dilating the mother wavelet (see Fig. [T|) , and 
convolving these daughter wavelets with the input sign al fsee Ens. [T2] and [T3|) . At a given location 
r and spatial scale s, the daughter wavelet is given by (jAddisonll2002l ) 



{r' 



\2l 



s 2s2 

and the daughter of the Fourier-transformed mother wavelet is given by 

{suj - Uo)^ 



■0r,s('^) 



(2s)VVV4exp 



-lujr 



(10) 



(11) 



3.2. The Wavelet Transform 



For an evenly sampled radial signal, x(r\ the continuous wavelet transform (CWtJ^ is gi 
by (jTorrence and CompolllQod : lAddisonll2002l ) 



iven 



/oo 
x{r')il)*g{r')dr' , 
-oo 



(12) 



where * denotes the complex conjugate. In practice, the wavelet transform is more efficiently 
calculated in terms of the product of the Fourier transforms, 



r(r,s) 



(13) 



For more discussion, see iFargd (j 19921 ) and iTorrence and Compol (jl998l ). the latter of whose 
softwar^ we use for this purpose. 



3.3. Choosing the Central Frequency 

Wavelet analysis follows a form of the Heisenberg uncertainty principle, in that an increase 
in spatial resolution comes only at the cost of decreased frequency resolution, and vice versa (see 



■''We use a discretized form of the CWT, given the digital nature of our signal. This is not to be confused with the discret e 
wavelet transform (DWT), which is based on fundamentally differen t mathematics ijPaubechieslllQoj ^lGraps lhQQStlMaUatillQQSl) . 
For analysis of a feature in Saturn's rings using DWT methods, see lBendiova et al. I 1I1993I ') and lSpah n et al. 

^Software available at |http: //atoc. Colorado. edu/research/wavelets/| 
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Fig. 1. — The Morlet wavelet with ujq = 6. The sohd hne shows the real part, the dashed line the 
imaginary part. For ease of viewing, the gaussian envelope within which the wavelet oscillates is 
shown as a dotted line. The wavelet translates and dilates in order to interrogate the signal at all 
scales (s) and locations (r). 



Addisonll2002l ) . Classical Fourier analysis can be seen as an end-member in this trade-off, providing 
the optimum possible frequency resolution, but at the cost of having no spatial resolution at all. 
For the Morlet wavelet, this trade-off is governed by the central frequency ujq in Eqs. [9] to [TTJ 
Higher values of cause the Gaussian envelope to contain more oscillations of the sinusoid, and 
consequently yield higher spectral resolution at the expense of lower spatial resolution (see Fig. [16] 
for an example). In this paper, unless otherwise specified, we use a value of wq = 6 (although there 
is nothing intrinsically special about any particular value of wq, which need not even be an integer). 
A typical Morlet wavelet can be seen in Fig. [H 

We should note that the expression for the Morlet wavelet given in Eqs. [9] to [11] is only valid 
for Wo ^ 5. For smaller central frequencies, one cannot ignore the addition of a real- valued constant 
to the sinusoid, which is neces sary to ensure compliance with the requirement that an admissible 



wavelet must have zero mean (jFarge 



1992 



Addisonll2002l ). 



3.4. Wavelength and Wavenumber 

The scale parameter, s, is directly proportional to the more familiar Fourie r wavelength, Xp. 
For the Morlet wavelet, this relationship is given by (iTorrence and Compdll998l ) 



Atts 



(14) 
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The wavelet plots in this paper are shown with a y-axis that increases linearly with the wavenumber 
k = li^jXp (and thus inversely with s), as the trace of a density wave appears linear in that 
representation (see Eq. [5]). 



3.5. Inverse Wavelet Transform 



-ds. 



The wavelet transform can be inverted, recovering the input signal to good accuracy. Before 
inverting, a powerful method of filtering can be implemented by setting undesired elements in the 
transform array to zero. The inverse wavelet transform can be written (|Fargdll992l ) 

"'^''^ ~ Cs Jo .3/2 

The normalization constant is 

^ Jo ^ 

where V'o,i('^) is given by Eq. [11] with r = and s = 1. If the original input signal x{r) was real, 
simply take the real part of Eq. [TSl 



-da;. 



(15) 



(16) 



The constant Cs is calculated only once for a given mother wavelet (i. e., i n the case of the 
Morlet, for a given value of ujq). The discretized form of Torrence and Compo ( 19981 ) includes a 
different method of calculating C^, which we will not restate here. Using their method, we find 
empirically for ojq > 6, 

0.776 (^) , (17) 



TC98 



6 



with the uncertain ty in the last decimal place. This result allows the inverse wavelet transform to 
be calculated with iTorrence and Compd 's software, even when ujq ^ 6. 



3.6. Wavelet Energy and Phase 

Because our wavelet transform is complex-valued, it can be expressed in terms of two 2-D 
arrays: the wavelet energy spectrum (or "scalogram") 

Ew{r,s) = \T{r,s)\^ (18) 

and the wavelet phase 

(^e{tf) ^ 

Most plots of wavelet transforms in this paper will show only the energy scalogram, which gives 
the relative importance of different frequencies at a given location. 

The two-dimensional array of wavelet phases can be more usefully reduced to the average 
wavelet phase, (j)w(r), weighted by the modulus, over all scales s at a single location r. For 
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a quasi-sinusoidal signal, (pwif) gives the local proximity to peaks and troughs, even when the 
sinusoid's wavelength is inconstant (as it is for density waves), and can be compared to a density 
wave's instantaneous phase (Eq.[6|). As with the filtering method described above, selected values of 
T(r, s) can be set to zero before performing this calculation, in order to isolate a desired component 
of the signal. 

Considering each complex element of T{r,s) as a vector in the real-imaginary plane, (pwii^) 
can be calculated by finding the average vector, then taking its phase. This is done by separately 
summing the real and imaginary parts of T(r, s) over the desired scales: 

keeping in mind that vr/2 < (pwi"!^) ^ 3tt/2 iff the real part is < 0. 



3.7. Wavelet Ridges 



It can be useful to plot local maxima in the 2-D representation of wavelet information. Wavelet 
ridges are those locations in the (r, s) plane for which 

d[Ew{r,s)] 



ds 



0, 



(21) 



and are us eful for pinpointing important scales (and thus frequencies, from Eq. [T^ at a given 
location r ([Addison et aLll2002l ). 



3.8. The Cone of Influence 



The cone of influence (COI) is the region in the wavelet transform within which edge effects 
may contaminate the signal. The size of the COI varies with the scale s, as a consequence of 
the dilation of the mother w avelet. It is up to the user to de fi ne a COI with which he or she 
is comfortable ([Addisonl |2002| ) . We follow iTorrence and Compd (|l998l ) in using the distance over 
which the wavelet energy contributed by the edge discontinuity decreases by a factor of e, which is 
\/2s for the Morlet wavelet. 

This form of the COI is, more generally, the interval over which the wavelet signal from any 
localized feature would be expected to spread. For example, the wavelet transforms of many density 
waves in this paper show some signal extending to the left of r^, a smearing effect which is the 
inevitable consequence of the spatial-frequency trade-offs discussed in Section 13.31 The reader can 
use the COI to visualize the expected spatial resolution along any horizontal line in a transform 
plot. 

For a given Fourier wavelength Xp, the COI width ^/2s is proportional to the denominator in 
Eq. [T4l Thus, the COI is broader for higher values of ujq (see Fig. [T6|) . 



Table 1: Parameter values (see Eq. [T]) for synthetic density wave discussed in Section HI 



Parameter Name 


Symbol 


Value 


Resonance label 


m 


10 


Background surface density 


o-Q 


40 g/cm^ 


Resonance location 


Tl 


130,000 km 


Initial phase 


't>0 





Damping parameter 




10 


Amplitude 


Al 


1 



In all wavelet plots in this paper, the COI is denoted by a region filled with cross-hatching. It 
is sometimes absent, in cases where the radial scan has been cropped so that true edges are not 
shown. 



3.9. Significance Levels 



When uncertainties in the input signal are characterized by gaussian white noise (with standard 
deviation ainput), t he normalized wavelet ene rgy Ewir, s)/af^p^^ is distributed as with two 
degrees of freedom (iTorrence and Compdll998l ). For our radial scans, each element is individually 
calculated as the weighted ave rage of a number of image pixels, and thus possesses a roughly 
gaussian error estimate (Tinput (jPaper 3)8 Therefore, in this paper, the lowest level shown in 
each contour plot corresponds to values of Ew{r,s) equal to 

^input' Further contour levels are 

logarithmic, with three levels per order of magnitude in Ew- 



4. Wavelet Analysis of a Synthetic Density Wave 

We here describe a multi-step semi-automated algorithm, using wavelet methods to measure 
the parameters describing a density wave. To illustrate the utility and limitations of our method, we 
consider a signal whose characteristics are known beforehand. We use the density wave curve taken 
from Eq. [H with parameters set equal to round numbers typical of the mid-A Ring (Table [T|). We 
note that these values allow the dimensionless radial parameter to be calculated as ^ = C(r — rj^). 



''For a radial scan element x, which is a weighted average of pi xel values Xj with weights Wi, the uncertainty can be estimated 
with a form of the standard error of the mean (s.e.m), given by llBevingtonlll969l ) 

9 _ - xf] 



The methodology of lPaper 11 does not allow for the above calculation; instead, we estimate ainput as 2 times the unweighted 
standard deviation, normalized by the square root of the local width of the scan area in pixels, which we find to be a good 
approximation for the scans discussed in this paper. Since the only current role for dinp^t is to normalize contour plots of the 
wavelet transform, we judge this approximation to be sufficient for our purposes. 
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r - (km) 

Fig. 2. — a) Synthetic density- wave radial profile generated by Eq. [T]with the values of Table [TJ b) 
Wavelet phase, c) Wavelet energy. As with all wavelet plots in this paper, contours are logarithmic, 
with three contours to an order of magnitude, and the lowest contour level is at the square of the 
median error estimate for the scan (here we artifically set a input = 0.1). The region filled with 
cross-hatching is the cone of influence. The dashed line is the foreknown wavenumber kowif) 
(Eq. [5|), while the solid line shows the calculated location of a wavelet ridge (Eq. [2T]) : here, the two 
diverge visibly only for r — r/, < 40 km, corresponding to ^ < 6. 

where C = 0.147 km^^ (Eq. [2]). 

4.1. Obtaining the Density Wave's High- Frequency Shape 

Fig. [2] shows the radial trace of the synthetic density wave, A(T(r), from Eq. [T]with the values 
of Table dl along with the phase and energy of its wavelet transform. Note that A(j(r) as discussed 
in this section is identical with the generic x{r) discussed in Section [3l The phase (i>w{f) in Fig- Wp 
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Fig. 3. — a) Calculated wavelet ridge (Section l3.7p . as in Fig. [2b (solid line); the expected wavenum- 
ber kowi''') (dotted line). The region of the wavelet ridge plotted in bold was used in a linear fit 
(dashed line); the y = point is plotted as a solid diamond, b) All three curves from Fig. [3^, 
shown as residuals with the expected wavenumber kjjwir). c) Calculated wavelet phase cpwir), as 
in Fig. [2)3, "unwrapped" to show how phase accumulates quadratically (solid line); the expected 
phase (j)Dwii^) (dotted line). The region of (pwi''') shown in bold was used in a quadratic fit (dashed 
line); the zero-derivative point is plotted as a solid diamond, d) All three curves from Fig. [3h, 
shown as residuals with the expected phase <?!>Diy(r'). e) The input synthetic density wave from 
Fig. [2^ (solid line); the fitted density wave, after the analysis of Section [4. 11 but still with randomly 
chosen values of and Al (bold solid line). Fitted values given in the figure can be compared 
with the input parameters (Table [I]) used to generated the wave. 
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is 0° at local maxima of Aa{r), and 180° at local minima, as we would expect. Fig. [2}: also shows 
the relationship between the foreknown wavenumber (Eq. [5]) and the wavelet ridge (Eq. [2T]) . For 
the first few wavecycles, the ridge is drawn away by non-sinusoidal components in the signal, but 
farther downstream it agrees with the expected wavenumber. The residual between the two is 
shown in Fig. [3)3. A linear fit to the portion of the wavelet ridge that is closest to linear can be 
used to re-obtain the background surface density o"o and the resonance location from Eq. [5j the 
fitted values are also given in Fig. [Sh, with error estimates taken from the linear regression. 

Fig. [3t shows how the wavelet phase (j)wif) can be "unwrapped" by adding 360° at every 
successive "wraparound", to reveal a shape very close to quadratic. The residual between the 
wavelet phase and the expected phase (j)Dw{''') (Eq. E]) is shown in Fig. [Sji. As with the wavelet 
ridge, the regions closest to ri do not correspond well to the foreknown values. But the phase 
reaches stable behavior much sooner (at much lower values of r — ri) than does the ridge. As in the 
previous case, the portion that is closest to quadratic can be fit to re-obtain the initial phase (po, 
the background surface density ctq and the resonance location r^] the fitted values are also given 
in Fig. [3li, with error estimates taken from the linear regression. 

We use the quadratic fit to the phase, rather than the linear fit to the wavenumber ridge, not 
only because it obtains all three parameters at once, but also because more of the curve is available 
for fitting. In practice, it is important to look at the residuals for each individual fit, and judge 
those portions of the phase curve deviating from quadratic behavior, in order to exclude them; only 
portions that conform to quadratic behavior should be used in the fit. 

When this phase of the analysis is completed, the high-frequency behavior of the wave has 
been obtained, but the shape and height of the wave's envelope is still unknown, as illustrated in 

Fig. 



4.2. Obtaining the Shape of the Density Wave's Envelope 

Now that we have determined the three parameters that define the shape of the wave's high- 
frequency components ((Tq, r^, and (po), we now accept that shape as given and proceed to fit the 
two remaining parameters (A^ and S^d) that define the wave envelope (i.e., the amplitude and decay 
of the high-frequency components). 

In the absence of other effects, the shape of the amplitude modulation (i. e., the location of 
the point of maximum amplitude) is governed by the damping parameter S^d. In fact, the non- 
oscillatory components of Eq. [T]are proportional to ^e~^^^^^^^ . If a reliable representation of the 
amplitude can be obtained, the best way to fit for is to find the dimensionless radial parameter 
(Eq. [2|) at which its derivative is zero. Denoting this point as ^max , we then have 

Cd = (23) 
Fig. |4^ illustrates the result of this calculation, in which we take the absolute value of the input 
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r - (km) 

Fig. 4. — a) The absolute value of A(j(r) (solid line); the same curve, smoothed three times with a 
22-km boxcar filter, giving the low-frequency shape of the wave envelope (bold solid line); maximum 
point of the latter (vertical bold dashed line), b) The input synthetic density wave from Fig. [2t^ 
(solid line); the fitted density wave, after the analysis of Section [4.21 but still with randomly chosen 
value Al (bold solid line), c) The input synthetic density wave from Fig. [2^ (solid line); the 
final fitted density wave, after the analysis of Section [4.31 Fitted values given in the figure can be 
compared with the input parameters (Table [1]) used to generated the wave. 



waveform Aa{r) (after first ensuring that the waveform oscillates about zero) and smooth it three 
times with a boxcar filter. The boxcar's width should be > 3 times the largest peak-to-peak 
wavelength in the scan. 

When this phase of the analysis is completed, only the overall amplitude of the wave remains 
unknown, as illustrated in Fig. [Da. 

The effects of are observable only in waves that can be described throughout by the linear 
theory of Eq. [H and indeed the analysis in this paper is limited to such waves. By contrast, 
the stronger resonances in Saturn's rings (including nearly all of the waves observed by Voyager) 
quickly become non-linear — that is, the oscillations Act become comparable to the background 



surface density ctq, causing the assumptions underlying Eq. [T]to break down (see lShu et aLlll985l : 
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Borderies et aLlll986l ). Non-linear density waves are characterized by sharp peaks and flat troughs, 
though the wavenumber knwif) is still fairly well-described by linear theory. 



Consequently, has sometimes been ignored in previous work. Firstly, iLissauer and Cuzzi 



(|l982l ) discuss and calculate the fractional distance Xi\il at which a density wave becomes non- 
linear, a calculation wh ich assumes t hat th e amplitude grows linearly with no modulation from 
the term. Secondly, Rosen et al. ( 1991a ) entirely ignored the e~^^^^°^^ term in Eq. [T]in their 
fits of density waves, which they justified by noting that they fit only the first few wavecycles of 
their waves. It is important to recognize that both of these treatments apply only to strongly 
non- linear waves, and are too simplistic for the myriad of weaker waves observable with CassinVs 
higher resolution. 



4.3. Obtaining the Density Wave's Amplitude 



Finally, we calculate Al using a Levenburg-Marquardt least-squares fit (jPress et adll992l ) to 
compare the density wave's radial trace to the shape defined by do, tl, cpQ, and obtained using 
the processes described above. Fig. shows the fitted curve overlying the original synthetic wave, 
with the fitted parameter value also given. 

We note that it is possible to obtain Al and S^o simultaneously with the least-squares fit, 
comparing the density wave's radial trace to the shape defined by only ctq, ri, and (pQ. Although 
this method yields acceptable results, our analysis of the synthetic wave shows that obtaining 
independently is superior. Still, the two-parameter least-squares fit is occasionally necessary, such 
as when the wave is truncated by an image edge or by unrelated nearby radial structure. 

When the amplitude Al is given in terms of the perturbation in surface density, the mass of the 
perturbing moon can be obtained through approp ri ate coeffi cients from the Disturbing Function 



( Goldreich and Tremaine 1982 : Rosen et al. 1991a jbl. IPaper I ). Unfortunately, conversion from the 



measured brightness (I/F) to optical depth, and thence to surface density, is greatly complicated 
by the presence of self-gravity wakes and by the ambiguities inherent in viewing the unlit side of 
the rings (see Section [5.80 . As the photometric calibration is still incomplete, this paper focuses on 
the results obtained from the spatial information in the data (we do assume that brightness varies 
relatively smoothly with surface density, such that we can correctly identify maxima and minima 
in the latter). We leave a more thorough analysis of wave amplitudes to a future paper. 



4.4. Fitting the Parameters of Real Density Waves 



We apply the method described in this section to the SOI imaging data presented in this paper 
and in lPaper j . calculating parameter values for many density waves in the rings. Two modifications 
of the method are required to allow for the differences between our synthetic wave and the real 
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data. 

Firstly, we remove low-frequency components of the signal by setting selected elements of 
r(r, s) to zero, taking the inverse wavelet transform, and then starting the analysis with this 
filtered wave. This ensures that waves oscillate about zero (as does Eq. [1]). In most cases, it is easy 
to draw the dividing line through regions of the wavelet transform with essentially zero energy (for 
examples, see Figs. [6l and [THj) . More problematic cases are noted in Table [3l We do not filter out 
high-frequency noise, which is difficult to separate quantitatively from the desired signal. 

In thus filtering the signal, we explicitly assume that variations in the background surface 
density are either absent or unimportant in determining the wave structure. In practice, we find 
that even waves with significant variation in background brightness can be well fit to Eq. [H thus 
justifying this assumption. Waves for which such brightness variations are significant are noted in 
Table [31 

Secondly, we find in practice that fitted values of the initial phase 0o are too uncertain to 
be useful. They turn out to be highly c orrelated w ith r^, as the distances between and the 



peaks of the wave depend linearly on (pQ (IShulll984l ). Furthermore, it is easy enough to calculate 



the value of (po from the location of the perturbing moon on its orbit (the uncertainty of which is 
exceedingly small). Therefore, we constrain (pQ to maintain the expected value while fitting for the 
other parameters. 



4.5. Error Estimates 

Meaningful uncertainties can be obtained for the parameters calculated in Section 14.11 — 
namely uo, r^, ar id (pn — by simply propagating the normalized errors obtained from the quadratic 



fit to 0iy(r) (see lPress et a/.lll992l ). However, this treatment assumes that the residuals will have 
an uncorrelated gaussian distribution about the fitted curve. In our case this is not true. Because 
the wavelet transform gives up some radial resolution in order to obtain frequency resolution (see 
Section l3.3p . the effective radial resolution of the wavelet phase (pwi''') is lower than that of the 
original scan. Thus, since we have used the same radial scale for all quantities, we have oversampled 
the wavelet phase so that adjacent pixels are highly correlated. 

We estimate the oversampling factor by examining the Fourier transform of the residual phase 
from Fig. [3ll. Indeed, as seen in Fig. [U very little power resides at scales within an order of 
magnitude of the nominal pixel scale. We define the "effective pixel scale" Ago as the scale above 
which 90% of the Fourier power resides. Thus, the calculated error estimates must be multiplied 
by a factor \/Xgo- This has been done for the synthetic error estimates quoted in Fig. [3l and will 
be done throughout this paper. 



We estimate the uncertainty in our fitted values of as being comparable to the largest 
peak-to-trough wavelength in the wave, since it is these wavecycles that are smoothed over in order 
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3 10 30 100 300 
Scale (pixels) 



Fig. 5. — a) The Fourier transform of the residual phase from Fig. [3]1. b) The fraction of Fourier 
power at scales A > A', for all A'. The dotted lines show that, for our synthetic wave, 90% of the 
Fourier power resides at scales greater than Ago = 22.7 pixels. 

to estimate the shape of the wave envelope (Section 14. 2p . This error estimate is converted from 
units of km to the dimensionless units of by inserting it in place of the quantity r — in Eq. [2j 

4.6. Method Summary 

The steps used in this paper to fit density waves in Saturn's rings are described above in detail, 
illustrated in Fig. [6l and briefly enumerated here: 

1. The wavelet transform T(r, s) is taken of the input signal Acj(r). 

2. A line is drawn through the 2-D wavelet transform, passing through local minima, separating 
the wave signature from lower-frequency background variations. 

3. Elements of the wavelet transform with wavenumber below the line are set to zero, yielding 
the filtered wavelet transform T'{r,s). 

4. The inverse wavelet transform of T'(r, s) is taken, yielding a high-pass-filtered signal Aa'{r). 

5. The average wavelet phase of T'{r,s) is taken, yielding (f)'y^{r). 
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Fig. 6. — The density wave fitting process, illustrated using the Pan 19:18 ILR density wave, 
observed in Cassini image N1467345975 (see Figs. [12] and [H]) . The two red vertical do tted line s 
indicate the interval used for the quadratic fit. From top to bottom: a) Radial scan (jPaper 11 ) . 
b) High-pass-filtered radial scan, with the final fitted wave shown in green, c) Wavelet transform 
of radial scan, with blue line indicating the filter boundary, and the green dashed line indicating 
the fitted wave's wavenumber. d) Unwrapped wavelet phase, with green dashed line indicating 
the quadratic fit and green open diamond the zero-derivative point, e) Residual wavelet phase, 
showing that the interval used for the fit is the interval in which the phase behaves quadratically. 
f) Wave amplitude, the local maximum of which (vertical dotted line) gives ^d; scale bar indicates 
the smoothing length of the boxcar filter. 
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6. A radial interval [ri,r2] is defined, and the phase (j)'y[r{r) on ri < r < r2 is fit to a quadratic 
function. The residual between (p'^yir) and the fit is inspected to verify that [ri,r2] well 
describes the interval that behaves quadratically. 

7. If the amplitude of Aa'{r) has a well-behaved shape (verified by visual inspection), then its 
maximum used to determine S^o- 

8. Taking the parameters fitted heretofore as given, a Levenburg-Marquardt least-squares fit of 
Aa'{r) is used to determine the remaining parameters (either Al only, or both Al and ^d)- 

Note that steps 2 and 6 explicitly require human input, while all steps are visually monitored to 
ensure that calculated curves are well-behaved. 



Observations 



In this section we de scribe sev eral highlights that emerge from wavelet analysis of the Cassini ISS 
radial scans of the rings (jPaper J ). The images from which the scans used in this paper are derived 
are listed in Table [2j A general overview of the rings can be found in Fig. [TJ showing the locations 
and context of the figures that follow. 



5.1. Pan density waves 



The Cassini ISS data set from S OI was the first to detect density waves excited by Pan, 
the moon embedded in the Encke Gap (jPorco et a/.ll2005l ): these waves are found in considerable 
quantity. Because of Pan's proximity, the azimuthal parameter m is limited only by the spacing 
between resonances becoming insufficient for a density wave to develop. Distinct density waves 
with m as large as 90 are observed. The central part of Fig. [8] (with radial trace in Fig. [9^) shows 
a region dominated by the Pan resonances 55:54 through 62:61; fitted wave models are also shown 
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Table 2: Observing information for images used in this paper. 
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" Measured from the direction of Saturn's north pole (ring- plane normal). Note that N1467351049 and 
N1467351539 view the lit face of the rings, the others the unlit face. 



in Fig. [9^. This region provides a cross-check on our error estimates, since adjacent waves can be 
assumed to have similar parameters. Fig. [9b plots the fitted parameters, showing that the error 
estimates well describe their mutual consistency. 

Pan density waves are also detected exterior to the Encke Gap, constituting the first outer 
Lindblad resonances (OLR) to be directly observed. There are fewer clear examples of these, as the 
outer A Ring becomes increasingly crowded with other density waves as orbital radius increases. 
Fig. [10] shows a region of the A Ring featuring Pan OLRs, which are unique among observed density 
waves in the rings in that they propagate inward. 



5.2. Atlas density waves 



The SOI data set also includes waves excited by Atlas (jPorco et a/.ll2005l ). These are fewer in 
number than Pan density waves, as Atlas' orbit is more distant, just beyond the edge of the A Ring. 
Given its small size . Atlas' mass is much more poorly constrained by orbital integrations (in fact. 



Spitale et al\ (|2006l ) are the first to constrain it at all), increasing the importance of measuring its 



mass from the amplitude of its density waves. 

Fig. [TT] shows the Atlas 5:4 density wave, which is located in a Cassini Division ringlet, and 



Fig. 8. — A portion of Cassini image N1467351539. A series of small density waves, mostly due 
to Pan, appear like grooves on a record between the Mimas 5:3 bending wave (lower left) and the 
Prometheus 13:12 and Mimas 5:3 density waves (upper right). See Fig. [9] for analysis. 



which gives the first direct measurement of surface density in that region (jPorco et a/.ll2005l ). Note 
that this wave is perched on top of a small increase in brightness, and that the wavetrain is truncated 
at the point where the brightness reaches its maximum. 



5.3. Second- and higher-order density waves 



Several second-order density waves excited by the larger ring moons Prometheus and Pandora 
are visible in our d ata. Unlike the first-order waves associa ted with these moons, most of which were 
noted by Voyager (IRosen et a/.lll991al : ISpilker et a?.ll2004l ) , the second-order waves are weak enough 
that their dispersion remains strictly linear (that is, the induced variations in surface density Aa 
are much smaller in magnitude than the background surface density ctq). Thus, these waves are 
much better approximated by Eq. [H making them much better suited for measurements of ring 
parameters. 



The mass of the perturbing moon can also be inferred from these weak waves, complementing 
the mass estimates that have been made using orbital integrations (jJacobson and FrenchI |200J; 
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Fig. 9. — a) Radial scan (jPaper J ) from Cassini image N1467351539 (see Fig. [8]), showing a string 
of density waves due to Pan ILRs 55:54 through 62:61. Fitted waves are shown in green, b) Three 
fitted parameters from these density waves, with the mean value plotted as a horizontal dotted 
line. The Pan 57:56 density wave cannot be simply fit because of interference from the nearby 
Janus 16:13 resonance. 



Renner et al\ l2005l : ISpitale et al\ 120061 ) , though this must await detailed photometric modeling 
that can absolutely relate observed brightness to surface density. 

Second-order and third-order waves excited by the co-orbital moons Janus and Epimetheus 
are also present in our data. Because of the interplay of perturbati ons from multiple moons, wi th 
variable orbits w hich periodically change the resonance locations ( Dermott and Murray 1981al lbl: 
Yoder et a/.lll983l ). these structures are generall y too comp lex to be fit using the methods described 
in this paper. They are discussed more fully in lPaper III . 
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Fig. 10. — Radial scan (jPaper J ) and wavelet transform from Cassini image N1467351049, including 
density waves due to Pan OLRs 58:59 through 73:74, which propagate inward because they are 
exterior to Pan. Density waves from other moons affect the signal at some locations. Dashed lines 
indicate model density wave traces, assuming a background surface density do = 15 g/cm^. 
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Fig. 11. — Radial scan (jPaper t ) and wavelet transform from Cassini image N1467345208, including 
the Atlas 5:4 density wave located in the Cassini Division. Dashed line indicates model density 
wave trace, using the fitted surface density a = 1.53 g/cm^. 



5.4. First-order density waves 

Our data show that all first-order density waves raised by the main ring-moons (i. e., excepting 
the small satellites Pan and Atlas) show quantitative signs of non-linearity (see Section \4.2\\ . The 
closest to linear are a handful of Prometheus and Pandora waves in the inner A Ring; in particular, 
Prometheus 9:8 (Figs. [12] and [T4]) appears to the naked eye as a beautiful textbook example of a 




iMiiiii 



Fig. 12. — A portion of Cassini image N1467345975, showing the nearly-archetypal Prometheus 9:8 
density wave. See Fig. [T3] for analysis. 



linear density wave (compare with the synthetic wave in Fig. [2^). However, upon closer examina- 
tion, even this wave has relatively flattened brightness maxima and relatively "peaky" brightness 
minima. In addition to being discernible in the signal, this "peakiness" is visible in the wavelet 
transform as a series of vertical streaks, as the signal becomes less like a modulated sinusoid and 
more like a series of ^-functions (the Fourier transform of which has equal power at all frequencies). 

Another indication of non-linearity is that the wavenumber deviates from a linear function of 
orbital radius, instead trending towards longer wavelengths very close to resonance, and shorter 
wavelengths downstream. This gives the wave a "concave up" appearance i n wavelet space, closely 



rese mbling the behavior of the waven umber in the non-linear wave models of lBorderies et ali (119861 ) 



and iLongaretti and BorderiesI (|l986l ) , who show that it is due to the wave itself modulating the 
background surface density. However, the middle part of the wave, especially for those waves that 
only barely enter the non-linear regime, may still reflect the background surface density. A few 
such waves, including Prometheus 9:8, are fit using the methods of this paper. We find that the 
residuals of the middle part of these waves, with respect to the quadratic fit, are low, and thus we 
list them as part of Table [3] — though, among the fitted values, only o"o and are likely to be very 
meaningful. 
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Fig. 13. — Radial scan (jPaper I ) and wavelet transform from Cassini image N1467346211. A 
seemingly chaotic region between the strong Prometheus 12:11 density wave (cyan dashed hne) 
and the very strong Mimas 5:3 bending wave (purple dashed line) is in fact populated by three 
distinct weak density waves: Pan 46:45 and 47:46, and Pandora 19:17. Dashed lines indicate model 
density wave traces, assuming a background surface density ctq = 37 g/cm^. 
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Fig. 14. — Radial scan (jPaper J ) and wavelet transform from Cassini image N1467345975 (see 
Fig. [T2]) . including the weak density waves Pan 19:18 and Pandora 15:13 perched on the inward 
slope of the stronger Prometheus 9:8. Dashed lines indicate model density wave traces, assuming 
a background surface density ctq = 35 g/cm^. 



Fig. 15. — A portion of Cassini image N1467346329, showing a "corduroy" pattern caused by 
moonlet wakes excited by Pan. Also seen are the Pandora 11:10 and Prometheus 15:14 density 
waves. See Fig. 1161 for analysis. 

5.5. Separating Waves in the Frequency Domain 

Fig. [I3] shows the region between the Prometheus 12:11 density wave and the Mimas 5:3 
bending wave. Although no coherent structure is apparent in the radial scan of this region, wavelet 
analysis clearly reveals three distinct smaller waves: two Pan ILRs and a second-order Pandora 
ILR. 

In Fig. [T^ wavelet analysis of the radial trace of the Prometheus 9:8 density wave reveals that 
its inward regions are convolved with the second-order Pandora 15:13 density wave, as well as the 
Pan 19:18 density wave. 

Although these waves are not well enough separated to be quantitatively fitted by the method 
described in Section SI they illustrate the ability of wavelet analysis to separate out and identify 
waves invisible to the naked eye. Yet another example is the identification of the very weak 
Janus 14:11 density wave in Fig. [16] 
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Fig. 16. — Radial scan ([Paper 11 ) and two wavelet transforms from Cassini image N1467346329 (see 
Fig. [T5|) . showing wakes excited by Pan along with several density waves. The upper wavelet plot, 
like all others in this paper, uses a central frequency ujq = 6; the lower plot uses ujq = 12, resulting 
in increased resolution in the spectral (y) dimension at the expense of smearing in the radial (x) 
dimension (see Section 13. 3p . Density waves are clearer in the upper plot, including the strong 
Pandora 11:10 and Prometheus 15:14 waves, but also the weak third-order Janus 17:14. Dashed 
lines indicate model density wave traces, assuming a background surface density ctq = 40 g/cm^. 
Moonlet wakes ex cited by Pan are clea rer in the lower plot. The three dotted lines denote the 
frequency profiles (jShowalter et al.lll986l ) of wakes that have traveled 342°, 702°, or 1062° in synodic 
longitude since their last encounters with Pan. 



5.6. Pan wakes 



The "corduroy" pattern pervading the ir nage in Fig. 1151 is caus ed by moonlet w a kes ex cited 
by Pan. Such wakes (described in detail by Showalter et al. ( 1986 ) and Horn et al. (jl996l )) are 
a fundamentally different process from waves in that they do not propagate, but arise from ring 
particle orbits organized coherently by encounters with nearby Pan. 

Wavelet analysis of the image is shown in Fig. [T6l with two different values of the central 
frequency ujq (Section 13. 3p bringing out different aspects of the structure. Like previous authors, 
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Fig. 17. — Radial scan (jPaper 11 ) and wavelet transform from Cassini image N1467344391, showing 
an unidentified wave- like feature in a non-plateau region of the C Ring. Although the Tethys 6:1 
ILR is likely too weak to excite observable perturbations, its location (fixed by a fiducial feature 
elsewhere in the image) is intriguing. Dashed line indicates model density wave trace, assuming a 
background surface density uo = 0.7 g/cm^ 



we note that wake wavelengths in the data tend to be a few percent longer than the model predicts, 
likely due to mutual gravitation among the ring particles. Secondly, not only can we verify the 
existence of a second-order wake (which is simply structure that has survived after more than one 
synodic period after its creation by Pan), but we also see signs of a third-order wake. Thus, wake 
structure is seen right down to wavelengths near the resolution limit, indicating that we still cannot 
place a limit on the remarkable resiliency of wake structures against damping. 



5.7. Unidentified Waves 



Several examples exist in the SOI data set of structures that appear to be waves but are 
unidentified. 

The first is shown in Fig. \T7\ taken from an image of a non-plateau region of the C Ring. The 
Tethys 6:1 ILR at 90,279 km presents itself as an intriguing explanation for this wave-like structure. 
The surface density obtained for an m = 2 wave is (Tq ~ 0.7 g/cm^, which is lower than any other 
direct measurement of ring surface density, but which may be consistent with the low optical depth 
of this region. Unlike most images considered in this paper, the ra dial scale here is con strained 
by a nearby circular feature at 90,405.7 km (designated "lER 24" bv lFrench et al\ (jl993l )). so the 
spatial correlation between the resonance and the observed feature is not artificial. However, the 
calculated strength (torque) of the Teth ys 6:1 IL R is several orders of magnitude lower than many 
resonances observed to excite no waves (jPaper S ). 
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Fig. 18. — Radial scan ([Paper j ) and wavelet transform from Cassini image N1467344627, showing 
an unidentified wave-like feature in the B Ring. Dashed line indicates model density wave trace, 
assuming a background surface density do = 60 g/cm^, a resonance location ri = 94,485 km, and 
a azimuthal parameter m = 4. Solid line indicates a high-pass filter boundary, and lower panel 
shows the radial scan obtained by inverting the filtered wavelet transform. 



Fig. [18] shows a second unidentified feature, at the inner edge of the relatively featurele ss 
region of the inner B Ring known as the "Flat Spot" (designated "lEG 83" by lFrench et ali (Il993l )). 
Superimposed upon the ~ 100 km undulation is a quasi-oscillatory signal whose frequency increases 
linearly with decreasing orbital radius, reminiscent of an inward-propagating density or bending 
wave. Assuming a background surface density (Tq ~ 60 g/cm^, obtained from the nearby Janus 2:1 
density wave, we estimate the azimuthal parameter m ~ 4. 



5.8. Photometric Regime 

Neither ring surface density nor optical depth is directly observed in imaging data; rather, these 
quantities inust be obtained through photometric modeling of the directly-observed brightness (e. g.. 



Pones et a/Jll993l ). This is in contrast to radio and stellar occultations, which directly measure 
the optical depth. Detailed models to obtain optical depths from Cassini imaging data are not yet 
complete, which is why this paper focuses on the spatial information in the data, leaving analysis of 
wave amplitudes for the future. For this purpose, all that is necessary is to know whether maxima 
(minima) in brightness correspond to maxima (minima) in surface density, or vice versa. The 
former is the "normal-contrast" regime; the latter is "reverse-contrast." 



Images of the illuminated face of the rings always have normal contrast. However, ambiguity 
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Fig. 19. — Ring brightness (I/F) as a function of normal optical depth (r), under the assumption 
of single scattering, using geometrical parameters appropriate for our highest-quality images (solar 
incidence angle 114.5°, emission angle 62.8°). Vertical dotted line indicates the location of T^nax' 



exists for the unilluminated face where, for example, no signal might mean either a totally clear 
region or a totally opaque region. This simple illustration indicates that the maximum brightness 
occurs at a finite optical depth Tmax, with normal contrast holding for optical depths r < Tmax but 
reverse contrast for r > Tmax- The value of Tmax can be estimated using the assumption of single 
scattering in a classical many-particle-thick ring; the dependence on optical depth is proportional 
to e""^/^' — e ""^/^, where n' and m are the cosines o f the solar incidence and the emission angles, 
respectively (jChandrasekhaii 119601 : ICuzzi et a/.l 11984 ). For the observing geometry of the highest- 
quality images considered here, Tmax ~ 0.43 (Fig. [T9|) . 



Mean backgrou nd optical depths for the A Ring were measured by Voyager at r > 0.5 (e.g. 
Esposito et a/.lll983l ). lea ding one to expect reverse contrast in th e unlit-side images of the A Ring. 
However, recent studies (jColwell et a/.l 120061 : iHedman et a/.ll2007l ) indicate that, rather than being 
a homogeneous sheet at the mean optical depth, the A Ring is instead an "intimate mixture" of 
nearly opaque self-gravity wakes with much more diffuse material (T^ap ~ 0.15) in the gaps between 
the wakes. The latter may well dominate the scattering of light into the camera, putting the images 
in the normal-contrast regime. A second argument in favor of the normal-contrast regime is that 
it results in lower error estimates in our model fitsU 

However, we adopt here the classical assumption of reverse contrast for this observing geom- 
etry. The morphology of non-linear density waves in the high-resolution images f Sectio n 15. 4D. the 



brightness minima of which are "peaky" while the brightness maxima are flattened (see IShu et al 



"The ratio in the uncertainties of ao and rj^ (the contrast ambiguity does not affect or Aj^) between reverse-contrast 
and normal-contrast ranges from 0.8 to 6.f , with a geometric mean of 1.5. 
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19851 ). is most naturally interpreted with the reverse -contrast assumpt ion. The morphology of 
the "straw" texture observed in strong density waves (|Porco et a/.ll2005l ) similarly favors reverse- 
contrast. Finally, the surface densities derived under the reverse-contrast assumption agree better 
with the work of previous authors, with waves from the lit-side (but lower-resolution) images re- 
ported here, and with preliminary analysis of more recent high-resolution lit-side data (not reported 
here) . 



6. Results of Wavelet Analysis 

The density wave parameters measured using the method described in Section S] are given in 
Table [3] and plotted in Figs. [20] and EH 



6.1. Surface Density 

Our measurements of the background surface density (Figs. [20k and l2ip are similar in magni- 
tude to previous measurements made using stronger waves in Voyager dat a (lEsposito et a/.[ [19831 : 
Lissauer[[l985[ : [Nicholson et a/.[[l99d : [Rosen et a/.l[l991al : [Spilker et a/.[[2004l ). but show significantly 
less scatter. We find a gentle linear trend of 0.0013 g/cm^/km through the inner and mid-A Ring 
(Fig. [2TI1 . while analysis of waves in the outer A Ring (which have too few wavecycles to be well 
characterized by our method) by other authors indicates that do decreases again exterior to the 
Encke Gap. The goodness-of-fit for this linear trend (reduced = H) with 29 degrees of freedom) 
suggests that the error estimates we quote for cjo are roughly 3 times too small, if indeed the linear 
trend is real. On the other hand, it may well be that at least some of the deviation from the linear 
trend reflects real structure in the A Ring. 

The weak waves (which are more likely to give well-behaved estimates of ring parameters 
- see Section 15. 3p in the unlit-side images (which have 3 to 4 times better resolution than the 
lit-side images - see Table [2|) hew most closely to the linear fit, with a standard deviation of 
2.1 g/cm^. Several of the significant outliers are relatively stronger waves, which may raise their 
own background surface density as very strong waves are known to do (Section [53]), indicating that 
the linear trend may be a more accurate reflection of the unperturbed surface density. 

A correlation can be observed between the size of the error bar and the damping parameter 
^D- This is because a lower value of denotes a smaller number of wavecycles before the wave 
damps away; consequently, the wavelength dispersion is less accurately measured, driving up the 
uncertainty in fio. 
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Table 3: Fitted values of density wave parameters including resonance location (r/,); deviation of the 
same from nominal pointing (Ar/,); background surface density (o"o); initial phase {(po); damping 
parameter (Cd)] wavelength, in pixels, at maximum amplitude (Xmax)', and amplitude (Al). Fits 
assume the reverse-contrast regime (Section 15. 8p . 





Image 


rL (km) 


ArL (km) 


O"0 


(g/cm^) 


<^o n 






\max (px) 


Al (I/F) 


Atlas 5:4 


N 1467345208 


118830, 


,62 


-9, 


,84±0, 


,22 


1, 


,52±0, 


.06 


323.1 


6, 


.97±1.47-'' 


8.0 


41.51 


Pan 7:6" 


N1467345326 


120669, 


,37 


-7, 


,65±0, 


,11 


3 


.32±0, 


.04 


201.2 


7, 


.55±0.85-'^ 


9.4 


15.68 


Atlas 7:6 


N1467345621 


124347, 


,51 


-4, 


,24±0, 


,30 


32, 


.61±0, 


.21 


0.4 


14, 


.99±0.60 


16.6 


9.10 


Pan 10:9'' 


N1467345621 


124609, 


,75 


-7, 


,45±0, 


,25 


37, 


.65±0, 


.42 


328.7 


11, 


.25±0.80 


19.5 


9.04 


Atlas 8:7^= 


N1467345739 


126048, 


,12 


-5, 


,61±0, 


,24 


35, 


.51±0, 


.21 


288.6 




13.38" 


18.9 


7.88 


Pan 12:11^^ 


N1467345739 


126126, 


,36 


-6, 


,09±0, 


,18 


36 


.70±0, 


.22 


293.3 




11.48" 


17.9 


9.80 


Pan 13:12 


N 1467345798 


126707, 


,39 


-6, 


,10±0, 


,11 


37, 


.51±0, 


.15 


274.8 


11, 


.60±0.72 


17.5 


7.59 


Prometheus 15: IS'^ 


N1467345798 


126800, 


,70 


-6, 


,72±0, 


,11 


38, 


.28±0, 


.17 


146.1 


11, 


.90±0.82 


16.6 


6.57 


Pandora 13:11 


N1467345798 


126897, 


,51 


-6, 


,26±0, 


,26 


37, 


.45±0, 


.34 


275.6 


12, 


.47±0.69 


17.1 


7.46 


Atlas 9:8^= 


N1467345857 


127363, 


,40 


-5, 


,30±0, 


,27 


37, 


.88±0, 


.27 


215.8 


10, 


.81±1.20 


23.8 


9.79 


Prometheus 8:7 


N1467345857 


127613, 


,59 


-10, 


,15±1, 


,14'* 


49, 


.08±1, 


.47 


178.8 


10, 


.26±1.15 


30.6 


173.90'' 


Atlas 10:9 


N1467345916 


128411, 


,03 


-6, 


,06±0, 


,35 


39 


.93±0, 


.51 


144.2 


9, 


.50±0.72 


26.9 


11.32 


Pan 19:18= 


N1467345975 


128897, 


,87 


-2, 


,10±0, 


,26 


39 


.43±0, 


.52 


167.6 


8, 


.33±1.03 


22.0 


13.63 


Prometheus 9:8 


N1467345975 


128945, 


,58 


-5, 


,02±0, 


,71'* 


46, 


.53±1, 


.02 


138.3 


9, 


.02±0.93 


33.1 


195.57'' 


Pan 22:21 


N1467346034 


129541, 


,69 


-2, 


,94±0, 


,30 


40, 


.87±0, 


.66 


113.8 


8, 


.14±1.25 


21.7 


16.28 


Pan 23:22 


N1467346034 


129718, 


,71 


-3, 


,83±0, 


,63 


41, 


.52±1, 


.87 


96.5 


8, 


.16±1.07 


21.4 


14.38 


Pandora 8:7 


N1467346034 


129746, 


,97 


-4, 


,70±0, 


,64'' 


41, 


.67±0, 


.70 


123.0 


9, 


.08±1.13 


34.2 


142.49'' 


Pan 24:23* 


N1467346034 


129880, 


,87 


-5, 


,80±0, 


,57 


47, 


.93±2, 


.23 


79.2 


7, 


.28±1.23 


25.3 


9.38 


Pan 27:26 


N1467346093 


130294, 


,86 


-4, 


,57±0, 


,44 


43, 


.49±1, 


.09 


24.6 


7, 


.74±1.30 


21.8 


12.67 


Pan 29:28 


N1467346093 


130523, 


,00 


-4, 


,82±0, 


,37 


40, 


.78±1, 


.14 


349.7 


7, 


.56±1.17 


20.9 


9.79 


Pan 34:33 


N1467346152 


130975, 


,36 


-2, 


,23±0, 


,47 


40, 


.16±1, 


.31 


259.1 


7, 


.60±0.80 


19.4 


9.64 


Pan 35:34 


N1467346152 


131050, 


,25 


-2, 


,71±0, 


,11 


40, 


.90±0, 


.36 


241.6 


7, 


.15±0.98 


20.5 


12.96 


Pandora 9:8 


N1467346152 


131101, 


,74 


-5, 


,83±1, 


,04'* 


42, 


.43±1, 


.42 


255.3 


8, 


.70±0.99 


35.5 


102.75'' 


Pan 42:41 


N1467351539 


131474, 


,25 


-4, 


,56±0, 


,90 


50, 


.58±4, 


.48 


118.4 


6, 


.30±1.25-'^ 


6.3 


9.04 


Prometheus 12:11 


N1467346211 


131590, 


,30 


-5, 


,68±0, 


,85'' 


48, 


.12±1, 


.71 


15.7 


7, 


.95±1.15 


36.1 


161.36'' 


Prometheus 12:11 


N1467351539 


131590, 


,30 


-3, 


,30±0, 


,92'' 


44, 


.10±1, 


.75 


324.2 


6, 


.96±1.10 


10.3 


147.91'' 


Pan 55:54 


N1467351539 


131974, 


,51 


-2, 


,91±0, 


,45 


40, 


.20±2, 


.05 


31.8 


5, 


.38±1.39-'^ 


5.8 


14.31 


Pan 56:55 


N1467351539 


132003, 


,34 


-3, 


,96±0, 


,37 


45, 


.77±1, 


.83 


52.8 


4, 


.58±1.14-'^ 


7.2 


25.55 


Pan 58:57 


N1467351539 


132058, 


,02 


-3, 


,40±0, 


,41 


43, 


.11±2, 


.12 


94.9 


4, 


.19±1.24-'^ 


7.5 


30.00 


Pan 59:58 


N1467351539 


132083, 


,96 


-3, 


,54±0, 


,54 


43, 


.90±3, 


.26 


115.9 


4, 


.71±1.39-'^ 


6.6 


24.41 


Pan 60:59 


N1467351539 


132109, 


,04 


-3, 


,33±0, 


,69 


42, 


.38±4, 


.67 


136.9 


4, 


.52±1.68-'^ 


6.7 


26.95 


Pan 61:60 


N1467351539 


132133, 


,29 


-2, 


,74±0, 


,38 


38, 


.99±2, 


.10 


158.0 


4, 


.46±1.61-'^ 


6.5 


22.42 


Pan 62:61 


N1467351539 


132156, 


,76 


-3, 


,31±0, 


,59 


43, 


.28±3, 


.85 


179.0 


3 


.62±1.48-'^ 


8.4 


45.14 



" The quality of the fit to this wave is affected by variations in background brightness (see Section I4.4|l . 
^ The quality of the fit to this wave is affected by the nearby edge of the image (see Section [4. 4p . 
'^ The quahty of the fit to this wave is affected by nearby density waves or other radial structure (see Section [4. 4p . 
'' This wave exhibits significant non-linearity ("peakiness"), which affects the quality of the fits to Atl and Al- 
For this wave, was fit simultaneously with Al, as described in Section [4.31 

At the location of this wave's maximum amplitude, the wavelength is less than 10 pixels, which may artificially 
lower the inferred (Section 16. 2|) . 
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Fig. 20. — a) Surface density ctq in Saturn's rings, as measured using density waves here and in 
previous work. Vertical lines indicate error bars; for points with no visible error bar, the error bar 
is smaller than the plotted symb ol. Error bars for pr evious work (not shown) are generally similar 
in magnitude to their scatter (see lSpilker et al.ll2004l ). b) Damping parameter ^/j in Saturn's rings. 
Vertical lines indicate error bars, c) Viscosity in Saturn's rings, as calculated here and in previous 
work. The associated a pparent ring th i cknes s is an upper limit (see text). Solid line indicates the 
theoretical prediction of lDaisaka et al\ (j200ll ). as calculated from Eq. | 
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Fig. 21. — Surface density ctq in the inner- and mid- A Ring for our results only (same data as in 
Fig. [20]) . Symbols indicate whether a wave appears in an unlit-side or a lit-side image (the former 
have better resolution), and whether an unlit-side wave is weak or strong (the former are more 
well-behaved). Dotted line shows linear fit ctq = (33.7 ± 0.1) + (1.3 ± 0.5)(riooo — 124), where cfq is 
in g/cm^ and riooo is in thousands of km. 



6.2. Damping Parameter and Viscosity 



Fig. 120b plots the damping parameter, giving a systematic survey of viscosity in the region 
between the Cassini Division and the Encke Gap. Fig. [20b plots the viscosity as derived from 
Eq. [71 showing a consistent increase over two orders of magnitude within the surveyed region . 



Again we have agreement with pre vious measurements (jEsposito et al\ Il983l : iLissauer et al 



Shu et al 



1984 



19851 : [Chakrabartil[l989l ). though the trend is much clearer in our data. We note that 



Esposito et ali s three innermost data points, which are discordant with our data, come from waves 
— Pan dora 5:4^ Prometheus 6:5, and Janus 4:3 — that are now known to have highly irregular 
shapes ([Paper j ). and thus to be unsuitable for such analysis. 



Caution must be exercised in measuring the damping parameter for a wave with wavelengths 
not much larger than the image resolution. Since a spiral density wave's wavelength decreases 
downstream (Section 12. 2p . the apparent amplitude can reach a local maximum because peaks and 
troughs are blurred as they become increasingly close together, not because of the actual physical 
damping process. Since is obtained by locating the point of maximum amplitude (Section l4.2p . 
the wavelength at that location can be found by combining Eqs. [H [3l [5l and [23] 



^1/3 ( 3("-^- 1) 
V27rG(Tori 



1/2 



(24) 



This characteristic wavelength (A^ 



max ) is given for each wave in Table [3] as a ratio with 
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the image resolution. Most waves in our data set are well-resolved, with generally 20 or more 
pixels per wavecycle at the maximum amplitude point. The waves in the lit-side images (Table [2]), 
however, are much more poorly resolved, with generally less than 10 pixels per wavecycle (keeping 
in mind also that the nominal image resolution is an upper limit; degradation can occur due to 
smearing, or imperfect pointing when taking radial scans). These waves also yield much higher 
viscosities than any previous author has reported for the A Ring. Given that poor resolution would 
tend to boost the apparent viscosity, we do not plot these waves in Fig. [20b and 120b. though we do 
report the values in Table El 

The rms velocity between ring particles can be derived from the viscosities using Eq. [8l and 
the vertical scale height of the rings is simply H ^ c/^ under the naive assumption that ra i idom 



veloc ities in the rings are isotropic. However, self-gravity wakes (IJulian and Toomre 



1966 



Sale 



19951 ) , which play a dominant role in driving viscosity in the A Ring (jPaisaka et al: '200 l|), cause 
rand om velocities to be larger within the radial plane than in the vertical direction (Daisaka and Ida 
19991 ). thus depressing the vertical scale height implied by a given magnitude of rms velocity. Hence, 
our "apparent thickness" should be interpreted as an upper limit for the A Ring. 

This "apparent ring thickness" can also be read from Fig. [20b . To allow viscosity and rms 
velocities to appear on the same plot, we used an average radial value to obtain fi, rather than 
using the actual orbital radius of each feature; this introduces errors of up to 4%. We also used an 
average optical depth r = 0.5. 

For the Cassini Division, however, where optical depths are too low for self-gravity wakes to 
be important, the implied vertical scale heights of 3.0 meters (Atlas 5:4) and 4.5 meters (Pan 7:6) 
are most likely correct as statedH The apparent thicknesses in the inner A Ring, consistently 
below 15 meters, also provide interesting constraints. However, the mid-A Ring's trend upwards 
to 35 meters near the Encke Gap probably gives larger-than-real thicknesses. 



Daisaka et al\ (|200ll ) parameterize ring viscosity from self-gravity wakes as observed in their 
simulations. Using Saturn's mass Ms = 5.69 x 10^^ kg and the density of ice p = 0.9 g/cm^, their 
expression reduces to 



v{r) ~ 26 



(25) 



122,000 km^ 

This estimate, which implicitly includes a prediction of the abundance of self-gravity wakes, is 
also plotted in Fig. [20b . Our data correspond fairly well, falling somewhat below the predicted 
viscosity in the inner A Ring (r < 127, 00 km), indicating that self-gravity wakes are somewhat 
less abundant there than [Daisaka et al\ predicted. The reverse is true in the mid-A Ring. We 
note that the peak of the a zimuthal brightnes s asymmetry, another indicator of self-gravity wakes, 
occurs at r ~ 130, 000 km ([Dones et a/.lll993l ). 



"As these tightly- wound waves have only ~ 10 pixels per wavecycle, their derived viscosities also may be anomalously high 
(see above). In this case, the implication is that the Cassini Division's vertical scale height could be even lower than reported 
here. 
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Table 4: Central pixel orbital radius for selected Cassini images, from fitted density wave parameters 
in Table [3l which assume the reverse-contrast regime (Section 15. 8p . For images containing more 
than one wave, the right-hand column shows the ratio between the uncertainty in rc and that 
predicted from the individual uncertainties in ri- 



Image 


rc (km)" 


^ waves 




N1467345208 


118735. 55±0.22 


1 




N1467345326 


120489.23±0.11 


1 




N1467345621 


124433.54±2.23 


2 


8.19 


N1467345739 


126012. 76±0.32 


2 


1.60 


N1467345798 


126749.02±0.37 


3 


2.84 


N1467345857 


127488.91±0.27 


1 




N1467345916 


128266.34±0.35 


1 




N1467345975 


129009.84±0.26 


1 




N1467346034 


129707.68±1.33 


3 


3.14 


N1467346093 


130398.98±0.18 


2 


0.44 


N1467346152 


131086.02±0.16 


2 


0.99 


N1467351539 


131976.81±0.52 


8 


1.08 



" Orbital radius associated with the central pixel [511.5,511.5] of the image 

6.3. Navigating Images 

Navigation data on the spacecraft's position and orientation, while excellent, are often accurate 
only to within tens of pixels for CassinVs Narrow Angle Camera (1 NAC pixel = 6 /.trad). Fiducial 
features appearing within an image are required for more accurate navigation, but truly circular 
features are rare in the A Ring. Here we use the density waves analyzed by our method as fiducial 
features to more accurately navigate the images in which they appear, by comparing the resonance 
locations determined by our method (i.e., where they are in the image) to the resonance locations 
that can be calculated from external parameters (i.e., where they should be in relation to Saturn 
and the rings). In images from the highest-resolution unlit images (Table [2]), lines of constant 
orbital radius appear nearly straight, so that fine navigation is only needed (or possible) in the 
radial direction. For these images, it is sufficient to identify the Saturn-centered orbital radius rc 
for the central pixel. This is done in Table HI 

Several images contain more than one density wave, providing a cross-check on the accuracy 
of our measurements. For such images, the quoted rc is the average weighted by the individual 
uncertainties, and the right-hand column in Table H] gives the re-normalization of the uncertainties 
as indicated by the agreement among the measurements. The average renormalization indicates 
that we have underestimated the uncertainty of rL by a factor of 3 (we came to a similar conclusion 
regarding our error estimates for ctq in Section 16. ip . 

Excluded from this analysis are all waves due to first-order resonances of the larger moons, 
whose deviations from linear theory (Section I5.4p affect the quality of the fits for r^. 
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7. Conclusions 

We have provided a primer on how wavelet analysis can be employed to study density waves 
and other radial structure in Saturn's rings. With this technique we have illuminated a number of 
ring features, as well as fitted the wave parameters of 32 density waves, most of them previously 
unobserved. 

Our results indicate a gently increasing trend of surface densities across the inner and mid- 
A Ring. We show viscous damping (and thus in-plane velocity dispersion) increases consistently 
from the Cassini Division to the Encke Gap, and place strict limits on the vertical thicknesses of 
the Cassini Division and the inner A Ring. We use the measured origin (resonance location) of each 
density wave to accurately navigate those images in which density waves appear that are suitable 
for our procedure. 

The final measurable parameter, the wave's amplitude, can be used to infer the mass of the 
perturbing moon, but only when the photometry is sufficiently well understood to relate observed 
brightness to ring optical depth and surface density. This will likely be carried out in a later paper. 

The 32 density waves that are quantitatively analyzed in this paper represent only a fraction 
of the previously unobserved density waves that can be discerned in Cassini ISS data using wavelet 
analysis. Some examples are shown here (see Figs. [TOl [131 fl^ and 1160 . in order to illustrate the 
power of wavelet analysis to identify waves that are very weak or packed closely together. We hope 
to undertake a more systematic survey of these in the future. 

We thank J. Harrington, C. Torrence and G. Compo for helpful discussions. We thank 
M. W. Evans for help with software development, and E. Baker for technical assistance with the 
radial scans. We thank J.-M. Petit, J. Cuzzi, and J. Scargle for significant improvement to the 
manuscript. We acknowledge funding by Cassini and by NASA PG&G. 
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